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Abstract. We apply the Bogomol'nyi technique, which is usually invoked in 
the study of solitons or models with topological invariants, to the case of elastic 
energy of vesicles. We show that spontaneous bending contribution caused by 
any deformation from metastable bending shapes falls in two distinct topological 
sets: shapes of spherical topology and shapes of non-spherical topology experience 
respectively a deviatoric bending contribution d la Fischer and a mean curvature 
bending contribution d la Helfrich. In other words, topology may be considered 
to describe bending phenomena. Besides, we calculate the bending energy per 
genus and the bending closure energy regardless of the shape of the vesicle. As 
an illustration we briefly consider geometrical frustration phenomena experienced 
by magnetically coated vesicles. 



PACS numbers: 02.40.-k, 87.16. Dg, 75.10.Hk, 11.27.-l-d 

Our motivation is to amplify on the observation of vesicles with arbitrary 
low genus (number of holes/handles) exhibiting conformal diffusion (spontaneous 
conformal transformation), namely existence of two conservation laws for vesicles 
[0-^. The genus is a topological invariant: a quantity conserved under smooth 
transformations which does not depend on the static or dynamic equations of the 
system under consideration. In contrast, conformal diffusion provides evidence that 
the system Hamiltonian is invariant under conformal transformations. Whereas the 
Noether theorem ^ || may be used to treat the latter invariance law in order to 
compute the corresponding conserved current and constant charge, the Bogomol'nyi 
technique enables one to treat successfully various models with topological invariants 
P-pd| . In this Letter we focus on the topological conservation law only; we defer the 
conformal diffusion to future articles. 

To obtain the Bogomol'nyi relationships we write down for vesicles of arbitrary 
genus a bending Hamiltonian as a covariant functional invariant under conformal 
transformations which depends on their shape only and which is suitable for the 
Bogomol'nyi decomposition. Applying the converse of the remarkable theorem of 
Gauss enables us to construct such an Hamiltonian. Instead of describing shapes 
by their Monge representation {i.e., their surface equation) as customary we 
characterize shapes by their metric tensor and their shape tensor: the total integral 
of the shape tensor self-product is our successful candidate. Then the Bogomol'nyi 



Bogomol 'nyi decomposition for vesicles of arbitrary genus 



2 



technique reveals the topological nature of bending phenomena while differential 
geometry extends forward the Bogomol'nyi relationships. In particular we show 
that any deformation of the non-trivial metastable shapes spontaneously leads to a 
deviatoric bending contribution a la Fischer p"^-p^ for shapes of spherical topology 
and to a mean curvature bending contribution (up to a conformal transformation of 
the ambient space) a la Helfrich for shapes of non-spherical topology: our 

approach exhibits that the bending contribution expression depends on the shape 



topology — our main result concisely contained in formulae (pO[), (21) and (p6|). 

Before describing the bending energy of vesicles, we first succinctly recall the 
Bogomol'nyi technique [Eqs. (0)-(|g) below] through the nonlinear cr-model |P,[l0|,[l7|]. 
More precisely we consider spin fields on curved surfaces S with the nonlinear cr-model 
as interaction |ll§-El|: 



^mag 



\jj^Ang'^Kpd,n''d,n^, (1) 
s 

where the order parameter n" corresponds to a point on the two-sphere §^ and the 
phenomenological parameter J to the coupling energy between neighbouring spins. 
The metric tensors gij and hap describe respectively the support manifold S {i.e., the 
underlying geometry) and the order parameter manifold §^ : as customary, g represents 
the determinant det{gij) and y^dfl the area element. Assuming homogeneous 
boundary conditions, when applicable, allows one to map each boundary to a point, 
then the support manifolds S are topologically equivalent to the torus Tg of genus g. 
Consequently, the order parameter field n" effects the mapping of the compactified 
support Tg to the two-sphere which is classified by the cohomotopy group 11^ (Tg) 
isomorphic to the set of integers Z : two spin configurations belonging to distinct 
classes cannot be smoothly deformed into one another. Since there exists an infinite 
number of classes that map the torus Tg of genus g to the two-sphere for each 
compactified surface S the space of spin configurations splits into an infinite number 
of distinct components, each characterized by a definite topological invariant. Here 
the topological invariant is the degree of mapping Q which is expressed in terms of 
the order parameter field as 

Q = 87 e*Va/3 (2) 
s 

with eij and fap the antisymmetric tensors associated with the support manifold 
S and the target (i.e., the order parameter) manifold S^, respectively. In spherical 
coordinates the total integral (^) immediately reads as the winding number. Note 
that the topological conservation law arises from the nature of the order parameter 
field only. Henceforth we demonstrate how the Bogomol'nyi technique enables one 
to study topological spin configurations subject to the magnetic Hamiltonian (|l|). By 
introducing the self-dual tensors 

T±,„ = ^ [9,rv T e„/„« S^-n^ , (3) 
which satisfy the precious inequalities P,p^ 

T±,„r±*" ^ 0, (4) 

the magnetic Hamiltonian (|l|) density decomposes as 

i9,n„5'n" = ir±,„T±^" ± fc.pd,n"d,nf . (5) 
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Integrating the previous decompositions (^) over the support manifold S and inserting 
the formula (||) lead us to rewrite the magnetic Hamiltonian (|^) in each topological 
class specified by the winding number Q in the form 

Timag = yJ^dQT^^aT^'" + inJ |Q|, (6) 

with (±) the sign of Q. Readily from the precious inequalities the magnetic 
Hamiltonian (|l|) yields 

nmag>4nJ\Q\. (7) 

Moreover, in each topological class, the lowest value of the magnetic energy is actually 
attained when the self-dual tensor T^ia vanishes. In other words, according to the 
definition (|^), the metastable order parameter fields n", which actually saturate the 
Bogomol'nyi bound (0), satisfy the first-order self-dual differential (or Bogomol'nyi) 
equation multiplet 

-±e;7"«<9,n^ (8) 

Clearly the Bogomol'nyi decomposition reveals the underlying topology: (i) the 
inequality (Q) claims that the minimum minimorum of the Hamiltonian in each 
topological class is proportional to the topological invariant; (ii) the metastable 
configurations which actually saturate the Hamiltonian satisfy something simpler 
than the usual (Euler-Lagrange) equations and no explicit solution is needed to 
compute their energy; (iii) from the computation emerges a self-dual symmetry; 
(iv) the decomposition clarifies the stability of the topological configurations and 
exhibits a spontaneous energy contribution due to any deformation from metastable 
configurations; (v) a deviatoric (or strain-like) tensor T^j^ is constructed. Usually 
not explored, the last two features will appear very relevant below. 

Now, we focus on vesicles of arbitrary genus. The large separation of length scales 
between the membrane thickness of vesicles and their overall size allows us to describe 
each vesicle as a surface manifold S embedded in the tridimensional space M?' [^,^. 
As a result, initially we focus only on the energy induced by the surface manifold S 
itself. On the other hand, the corresponding Hamiltonian has to be invariant under 
conformal transformations and suitable for the Bogomol'nyi technique. Nonetheless, 
we must first choose how to represent the surface manifold S: we invoke the converse of 
the remarkable theorem of Gauss Let S{x^) be a surface manifold with arbitrary 
intrinsic coordinates [x^) embedded in the tridimensional ambient space A'^, then 
there uniquely exist a metric tensor gij{x^) and a shape tensor hij{x^) related by a 
definite equation multiplet; the converse is true except for the exact position of the 
surface S in the ambient space A'^. Then, with a different notation, we concisely 
express the remarkable theorem as 

5(a:^)cA3 ^ (g,,(x'=),6,,(x'=)). (9) 

The geometrical meaning of the couple (g^ , bij) in (^) is revealed by introducing both 
the infinitesimal tangential displacement dT over the surface manifold S and the 
infinitesimal normal displacement dN to the surface manifold S inside the ambient 
space A^; we have 

dT-AT ^gij{x'')dx'dx^ , (10a) 
d^-dT = bij{x^)dx'dxK (10&) 



Bogomol 'nyi decomposition for vesicles of arbitrary genus 



4 



In differential geometry, the quadratic differential forms ( |lOq ) and (f06) are referred 
to as the first and second fundamental forms, respectively. Furthermore, the Riemann 
tensor, which measures how much a manifold is curved, reduces to 

Rklmn — K [gkm9in — 9kngim] , (11) 

where the Gaussian curvature K depends only on the metric tensor {gij) and its first 
and second derivatives: such an entity is called a bending invariant or is said to be 
intrinsic . When the surface manifold S is embedded in a curved tridimensional 

ambient space A'^, then the intrinsic curvature K splits as 

K^G + k, (12) 

where the extrinsic curvature G satisfies 

G=\e'^e^'hkhju (13) 

whereas the sectional curvature K characterizes the ambient space. For example, 
the sectional curvature K vanishes for the flat space M?" and is equal to one for 
the three-sphere S'^. We complete this brief overview by stating the Gauss-Bonnet 
theorem [|2^,^: 

j V5dr!/^ = -47r (g + e-1) (14) 

with iSg^e a surface manifold topologically equivalent to a closed surface manifold of 
genus g less e points (ends) embedded in A^^. 

With this background, we now assume the bending Hamiltonian as the following 
functional 

Hfc = ifc j ^dng'^gkib.H/, (15) 



where the phenomenological parameter k describes the bending rigidity. Since the total 
integral in (^) is known to be invariant under conformal transformations |25,^, it 
remains to show that our assumption fits the Bogomol'nyi technique as desired. To 
begin, note that the bending Hamiltonian Hb formula (jl^) stresses the similitude with 
the magnetic Hamiltonian Hmag (0)- Then the self-dual deviatoric tensors are defined 
as 

which yield the precious relationships 

B^^JB^'' = (^TXf ^0 (17) 

where A and A denote the eigenvalues of the shape tensor {b^j), namely the principal 
curvatures of the surface (or support) manifold S. Here the decompositions (||) read 

i6.,6''^- = iB%S±^^±G. (18) 

By inserting the previous decomposition ([l^) in the formula (|l^) and recognizing the 
total curvature (p^), the bending Hamiltonian Tib ( p^ ) decomposes as 

nb = ^k J V5df^S±y.S±'^ T47rfc (g+e-1). (19) 
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Note that we have used the fact that, according to (^_2|), the extrinsic curvature G and 
the intrinsic curvature K are equal when the ambient space is flat {K = 0). Existence 
theorems select the correct value for the sign (±) and the proper topological classes 
for vesicles as follows. From the precious relationships it is evident that the 
surface manifolds S which actually saturate the Bogomol'nyi decomposition (|l^) are 
the totally umbilical surfaces (A = A) and the minimal surfaces (A + A = 0). Since 
only the round two-sphere §^ is totally umbilical ||2^,|2^, the decomposition with sign 
(+) is relevant only for surface manifolds So topologically equivalent to the round 
two-sphere S^; our first key result reads 

Hb [5o] = 5^ / \/5df^ B+^^B+'' + 4Trk. (20) 

<SoCR3 

It is noticeable that any deformation of the metastable bending configurations for 
shapes of spherical topology spontaneously leads to a deviatoric bending contribution 
a la Fischer p^-p^: to the best of our knowledge, there is no direct derivation in the 
literature for such a bending contribution suggested first by Fischer . Since within 
the flat ambient space there is no closed minimal surface (e = 0) whereas minimal 
surfaces of genus g (g ^ 0) with e ends (e ^ 2) do exist |Q, the decomposition with 
sign (— ) is relevant only for surface manifolds Sg^e topologically equivalent to such 
minimal surfaces ^; our second key result immediately reads 

[5g,e] = iA; J ^dQB-^^B^'' +Ank{g+e-l). (21) 

Here any deformation of the metastable bending conflgurations spontaneously leads 
to a mean curvature (or anti-deviatoric) bending contribution a la Helfrich widely 
employed to describe bending phenomena Note that the minimal surface S_q 2, 

namely the catenoid, is the elementary neck used to build or analyze numerically 
vesicles of arbitrary genus g (^-^ . We cannot proceed much further unless we invoke 
an assertion due to Chen which claims ||25| , [2^ 

J^dfl^b,jV^= J ^dfl'\^^b',,b'''+K'Y (22) 

c(5)Cc(R3) 

where c corresponds to an arbitrary conformal transformation. Consequently the 
bending Hamiltonian Tib (^Sh extends as 



Tib^\k / ^'dn' b'.jb''^ +2K' . (23) 



c(5)Cc(R3) 

It is easily seen that the extended bending Hamiltonian Tib ( p3| ) decomposes as 
nb=\k J ^dn'B'^^^B'^"' T47^kig+e-l) + 2kCf 

with C+ = and = J y/g^ dn' K' . (24) 

c(5g.e)Cc(R3) 

Consequently, the Chen assertion (|2|) strongly suggests to apply existence theorems 
for minimal surfaces living in curved tridimensional space A'^. To illustrate this point, 
we will establish Bogomol'nyi relationships for closed surfaces. A theorem due to 
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Lawson claims that there exist closed minimal surfaces £,m,n of arbitrary genus 
g = mn in the three-sphere (K = 1); ^o,n is the round two-sphere ^i^i the flat 
torus {i.e., the Clifford torus). Furthermore, the Willmore-Kusner conjecture p8|-^ 
asserts that the surfaces ^g,i actually minimize the functional in (|2^): we may have 



Ag = inf 



c(5j,o)Cc(K3) 



V?df7. 



(25) 



Readily the Bogomol'nyi decomposition for closed surfaces (e — 0) takes the form 

rLb[S^]^\k J y^dn'B'-^^B'-'^ +ATTk{g-l) + 2kA^, (26) 

c(5g)CS3 

where the conformal transformation c maps M"^ to S"^. Accordingly, the Bogomol'nyi 
technique combined with the Chen assertion (|2|) allows us to measure how much 
a surface manifold is deformed from any surface which is minimal inside a certain 
curved ambient space A^. Thus, without loss of generality, the Willmore-Kusner 
conjecture ( |2^ ) enables us to outline through the formula ( p6| ) our third key result: 
any deformation of the metastable bending configurations for shapes of non-spherical 
topology spontaneously leads to a mean curvature bending contribution a la Helfrich 
up to a conformal transformation of the ambient space. Before summarizing let us 
expose how the bending energy bounds (i.e., the Bogomol'nyi bounds) are governed. 
From the Bogomol'nyi decompositions (|o|), (^ and ( p6| ) the Bogomol'nyi bounds 
are: 



'^&['5g,e] ^ 



2Agfe-H47rA;(g-l) 
Airk (g+e-1) 



if g^O and e = 0, 
if g^O and 8^2. 



(27) 
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Obviously the Bogomornyi bounds for vesicles ( p7| ) are proportional to the bending 
rigidity k and separate into two parts: a closure bending energy 2Agfc and a 
genus/end bending energy Ank (g+e— 1). Whereas the bending energy per genus/end 
is straightforward to compute {Ank), there exists no literal formula for the Lawson 
sequence Ag yet: for the round two-sphere and the flat torus the Lawson 
sequence Ag takes respectively as exact value Ao = 47r and Ai = 27r^ {circles in Fig. |l|), 
for higher genus we have computed numerical estimates {crosses in Fig. \]\). 

To summarize our key results concisely expressed in formulae ( pO| ) , and ( p6| ) 
the Bogomol'nyi technique, extended forward with existence theorems from differential 
geometry, allows us to show that any deformation of the non-trivial metastable shapes 
(which are clearly identified) generates a bending elastic energy contribution which 
falls in two distinct categories with respect to the topology of the shape: (i) for 
shapes of spherical topology a deviatoric bending contribution a, la Fischer |l3| |l5|] 
spontaneously arises; (ii) for shapes of non-spherical topology a mean curvature 
bending contribution a la Helfrich ||, ^|l6| (up to a conformal transformation of the 
ambient space) spontaneously emerges. We observe that this splitting comes from 
the breaking of the self-dual symmetry of the Bogomol'nyi relationships by existence 
theorems from differential geometry. Besides, the Bogomol'nyi relationships allow us 
to compute both the bending energy per genus/end and the bending closure energy 
for vesicles regardless of their shape. 

As a further illustration, our approach leads to a clear understanding of 
geometrical frustration phenomena experienced by magnetically coated vesicles [^8|- 
^ : the presence of the double Bogomol'nyi decomposition generates a competition 
between magnetic solitons and shapes which tend to saturate the magnetic energy (^ 
and the bending energy (p^), respectively. More precisely, when at least one deviatoric 
tensor cannot vanish, the balance between the two deviatoric energies releases the 
frustration: hence both magnetic and geometric effects manifest in accordance with 
the deviatoric energies, e.g. by removing a mismatch between magnetic and geometric 
(or underlying support) length scales. 

In conclusion, our approach gives in a rather natural manner the Bogomol'nyi 
relationships for vesicles: their typical features combined with existence theorems from 
differential geometry show that spontaneous bending deformation from metastable 
bending shapes splits in two distinct topological classes (shapes of spherical topology 
and shapes of non-spherical topology): in other words, topology may be considered 
to describe bending phenomena-in contradiction with customary phenomcnological 
approaches 1^, jsj . Furthermore, the appearance of Bogomol'nyi relationships for 
vesicles provides a powerful guide for understanding vesicles and enlarges the 
application field of the Bogomol'nyi technique (traditionally used in fields ranging 
from condensed matter physics to high energy physics) to elastic and geometrical 
phenomena in soft condensed matter. 

This work was supported in part by the U.S. Department of Energy. 
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